arXiv:1505.01468vl [cond-mat.other] 6 May 2015 


Transport in ultradilute solutions of ^He in superfluid '^He 


Gordon Baym,“’^ D. H. Beck,“ and C. J. Pethick“’^’° 

’^Department of Physics, University of Illinois, 1110 W. Green Street, Urbana, IL 61801 
^ The Niels Bohr International Academy, The Niels Bohr Institute, University of Copenhagen, 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
”NORDITA, KTH Royal Institute of Technology and Stockholm University, 
Roslagstullsbacken 23, SE-10691 Stoekholm, Sweden 
(Dated: May 7, 2015) 

We calculate the effect of a heat current on transporting ^He dissolved in superfluid '^He at ultralow 
concentration, as will be utilized in a proposed experimental search for the electric dipole moment of 
the neutron (nEDM). In this experiment, a phonon wind will generated to drive (partly depolarized) 
®He down a long pipe. In the regime of ®He concentrations and temperatures ~ 0.5 K, the 

phonons comprising the heat current are kept in a flowing local equilibrium by small angle phonon- 
phonon scattering, while they transfer momentum to the walls via the ^He first viscosity. On the 
other hand, the phonon wind drives the ^He out of local equilibrium via phonon-^He scattering. 
For temperatures below 0.5 K, both the phonon and ^He mean free paths can reach the centimeter 
scale, and we calculate the effects on the transport coefficients. We derive the relevant transport 
coefficients, the phonon thermal conductivity and the ^He diffusion constants from the Boltzmann 
equation. We calculate the effect of scattering from the walls of the pipe and show that it may 
be characterized by the average distance from points inside the pipe to the walls. The temporal 
evolution of the spatial distribution of the ^He atoms is determined by the time dependent ®He 
diffusion equation, which describes the competition between advection by the phonon wind and ®He 
diffusion. As a consequence of the thermal diffusivity being small compared with the ®He diffusivity, 
the scale height of the final ®He distribution is much smaller than that of the temperature gradient. 
We present exact solutions of the time dependent temperature and ®He distributions in terms of a 
complete set of normal modes. 

PACS numbers: 67.60.G- 13.40.Em 05.20.Dd 


I. INTRODUCTION 

The physics underlying the transport properties of 
mixtures of ^He and superfluid ^He changes markedly 
as the concentration of ^He varies. We determine here 
the transport properties of these mixtures at very low 
concentrations, X3 = n^/{n^ + n4),$10“®, where ng and 
714 are the ^He and ^He densities, and low temperatures, 
T fZjQ.Q K, where phonons are the dominant superfluid 
excitation. In this case, the phonons are in local ther¬ 
mal equilibrium; their interactions with the ^He distort 
the ^He distribution and dominate the ^He diffusion. For 
concentrations 2:3 ^ 10“^, the reverse situation holds: the 
^He are in local equilibrium due to rapid ^He-^He scat¬ 
tering and the phonon distribution is distorted due to 
phonon-^He interactions [I[. In the intermediate con¬ 
centration regime, the phonon and ^He distributions are 
both distorted and must be determined by solving the 
coupled evolution, or Boltzmann, equations 0 . At the 
highest concentrations, ~ 1 %, Fermi-Dirac statistics 
for the ^He become important Q. These transport prop¬ 
erties are of interest as an example of a two-component 
fluid with excitations of comparable energy but very dif¬ 
ferent momenta, and where the excitations of the two 
species obey different statistics. 

The transport properties of ^He in superfluid '^He at 
low concentrations are also important for the proposed 
experiment Q to measure the neutron electric dipole mo¬ 
ment (nEDM) at the Oak Ridge National Laboratory 


Spallation Neutron Source. There, the neutron preces¬ 
sion frequency will be determined using the absorption 
of polarized ultracold neutrons on polarized ^He atoms 
in solution in superffuid ^He via the reaction 

n -b ^He —>■ p -b t -b 764 keV, ( 1 ) 

which has a strong spin dependence, since capture pro¬ 
ceeds primarily through the spin-singlet channel. Two 
key considerations accrue from this choice of detection 
technique. In order to maximize the precision with which 
the precession frequency can be measured, the optimal 
^He concentration, x^ ^ 10“^°, corresponds to a cap¬ 
ture rate comparable to the decay rate of the neutrons. 
However, primarily due to wall collisions, the ^He will 
gradually become depolarized. In order to reduce the 
background from neutron capture on unpolarized ^He, it 
is crucial to be able to periodically sweep out the ^He by 
means of a heat current Q. In this paper, we calculate 
both the heat and ^He particle currents based on well- 
established microscopic theory of phonon-phonon Q and 
phonon-^He scatterings Q, as well as the evolution of 
both the temperature and ^He concentration. 

At the concentrations and temperatures of interest in 
the experiment, in addition to phonon-phonon, phonon- 
^He and ^He-^He scattering, the scattering of both 
phonons and ^He from the walls of the containers can 
also be important. Here we extend the solution of the 
Boltzmann equations in Ref. 0 to include these effects, 
in addition to providing some examples for xalSlO”®. For 
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illustration, we consider the effect of a heat current in an 
essentially one dimensional geometry, with the ^He su¬ 
perfluid ^He mixture in a long pipe with a diameter of 
a few cm. The phonon-wall interactions affect the ther¬ 
mal conductivity as well as the phonon velocity distribu¬ 
tion within the pipe; the ^He-wall interactions affect the 
transport of the ^He in the presence of a heat current. 

This paper is arranged as follows: Section II describes 
the basic scattering mechanisms the calculation of trans¬ 
port coefficients from the Boltzmann equation is given 
in Sec. III. Subsequently, we calculate the temporal and 
spatial evolution of the temperature (Sec. IV) and the 
^He density (Sec. V). We summarize results in Sec. VI. 
In Appendix A we analyze the transport when scattering 
of phonons is predominantly from the walls of the pipe, 
and in Appendix B we solve analytically the equation for 
the temporal evolution of the ^He concentration. 


II. PHONON AND ^HE RELAXATION 


We begin by considering the relevant microscopic re¬ 
laxation mechanisms (detailed in Ref. 0 ), first for the 
phonons. The momentum-dependent mean free path of 
a phonon of momentum q scattering against the ^He, 


^phz{q) 


s 

Iq 


47rn4 1 
x^J 


( 2 ) 


is typically greater than 1 km for xz rS 10“® and T ~ 
0.5 K 0; here 7 ^ is the corresponding scattering rate, s 
is the phonon velocity, and J is an angle-integrated rate 
constant. Therefore, phonon-wall and phonon-phonon 
scatterings determine the phonon distribution. As dis¬ 
cussed in Refs. [1] and [3, rapid, small angle phonon- 
phonon scattering establishes thermal equilibrium along 
phonon ‘rays,’ i.e., given directions in momentum space, 
with the distribution 


^9 ^{sq-q-Vph)/T{r) _ 

where T(f) is the local temperature and Vph is the mean 
phonon drift velocity. Large angle phonon-phonon scat¬ 
tering, either in a single event or a succession of small- 
angle processes, is slower and gives rise to the phonon 
first viscosity, 


riph = 7 -4 

5 s 


( 4 ) 


where Sph is the phonon entropy density and £msc is the 
viscous mean free path. At a pressure of 0.1 bar. 


£visc — cm, (5) 

to a good approximation i[i3 , where Tr- is the tem¬ 
perature in Kelvin; at T= 0.45 K, Ivisc — 0.17 cm. 


In the presence of a heat flux, Q = TSphVph, small- 
angle phonon-phonon scattering keeps the phonons in lo¬ 
cal thermal equilibrium, where the mean phonon drift 
velocity, Vph-, is 


TSphVph = -KphS/T. ( 6 ) 


The thermal conductivity of the phonons, Kph^ can, at 
low concentrations, be written as 


£^ph 


5 sSphR^ 

8 4 // 


( 7 ) 


in a pipe of radius R and where (.ef f is the effective mean 
free path 


^visc OIL 

The second term represents scattering of the phonons on 
the walls; the numerical coefficient 16/5 is chosen to give 
the correct Casimir limit. In this limit, ivisc large com¬ 
pared to the pipe diameter (see Appendix A), the phonon 
thermal conductivity assumes the Casimir form [ll| 

Rph,Casimir — ‘^RsSph- (9) 

In the opposite limit, ivisc ^ R, the thermal conductivity 
becomes 


K.. 


ph,visc 


5 sSphR^ 


( 10 ) 


The ^He contribution to the overall heat flux is negligible 
at low X 3 Hi- 

The mean free path of a ^He scattering on unpolarized 
^He is [1 


C33 = 


( 713 / 2 ) 0-33 


8.66 X 10 
X 3 


-8 


-cm. 


( 11 ) 


where 0-33 is the corresponding cross section. Thus for 
xz < 10“®, one has £33 4 1 m; we therefore neglect ^He- 
^He scattering. On the other hand, the mean free path 
for ^He scattering on phonons 0 , 


i^ph 


•\/3 / 774 \ ^ 77l*^/^S^ 
T3/2 


0.077 


0.45 

—— cm. 
T ) 


( 12 ) 


where m* is the ^He effective mass in superfluid ‘^He, is 
small compared to the pipe diameter for 0.3 K. Thus 
the dominant process for bringing the ^He toward equi¬ 
librium for temperatures of interest in the experiment 
is scattering against phonons. In the next section, we 
outline the calculation of the ^He transport coefficients; 
more details can be found in Ref. 0. 
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III. ®He BOLTZMANN EQUATION AND 
TRANSPORT COEFFICIENTS 


The ^He Boltzmann equation has the general form 

dfp , P 


.v/ 

ot m* 


p',q,q' 

-fpn^S{r){l + n^S,{r)] 
Sfp - /3fpV3 - f3f°V3 Sfp 


T 33 


T3u 


'^3wd 


Here fp is the ^He distribution function, 


(13) 


jO _ g-/3(p^/2m*-^i3) 


(14) 


is the equilibrium distribution function, and we write the 
deviations from local equilibrium as 

Sfp = fp-fp^°, 

where 

is the local equilibrium distribution function, i.e., the dis¬ 
tribution towards which collisions with phonons drive the 
^He. The first term on the right represents the scatter¬ 
ing on the phonons, the second term the scattering from 
other ^He (numerically insignificant for the concentra¬ 
tions of interest) and the last two terms isotropic dif¬ 
fuse (d) and specular (s) scattering of the ^He from the 
walls. In the term describing collisions with phonons, 
p and p' are the initial and final ^He momenta, re¬ 
spectively, and q and q' are the corresponding phonon 
momenta. The phonon~^He scattering kernel is T = 
\{p'q'\T\pq)\'^2TTS{p^/2m* + sq — p''^/2m* — sq'), and mo¬ 
mentum conservation, p' + q' = p + q, is understood in 
the collision term. 

In order to calculate the effects of a phonon wind on 
the ^He we solve Eq. (ICT for Sfp. On the left side of 
the Boltzmann equation, we approximate the distribu¬ 
tion by its local equilibrium form, We neglect the 

contribution from dvph /dz because of the relatively small 
temperature gradient (see Eq. (l38ll below), while the gra¬ 
dient of P in this term gives a second order contribution, 
which we neglect. The left side of the Boltzmann equa¬ 
tion is then 


(15) 

(16) 


dz 


P 


- -T 
2 m* 2 


/ 


.Q 1 ,.n 1 dn. 


pt 2 dz ■'P' 


+ fp-^- (17) 


On the right side we write 

Sfp= dfpPzWp] (18) 


To solve for Wp we multiply the Boltzmann equation by 
Pz and integrate over all p. Noting that the distortion Sfp 
is proportional to Pz, we see first that on the right side 
of the Boltzmann equation the two term oc Sfp — fifpV 3 
do not contribute, since both the ^He-^He scattering and 
the specular ^He scattering from the walls conserve mo¬ 
mentum in the z direction. Following Eq. (81) of Ref. 0 
for the phonon-^He scattering, we see that the remaining 
terms on the right side comprise 

- ^ - PfpPzVph] - (19) 

^3wd 

With the inclusion of the recoil effect in the phonon-^He 
scattering to lowest order (see the Appendix of Ref. [l| 
and Sec. 6 of Ref. i) , the solution of the Boltzmann 
equation is 


5fp 




fpVph 


df^ 

dz 


( 20 ) 


where 

1 _ d^rec ^ 1 

T 3 ~ m* 


( 21 ) 


is the effective ^He scattering rate, including both scat¬ 
tering from phonons, encoded in Tree, and diffuse scat¬ 
tering from the walls of the pipe. 

Integrating the Boltzmann equation, Eq. m, over p 
we recover the continuity equation 


dn 3 

^ + V.,3 = 0, 


( 22 ) 


with the ^He particle current given by 


js = V 


/ 


d^p p 
(27r) m* 


(23) 


where v is the number of spin degrees of freedom of the 
^He: 1 for a fully polarized sample and 2 in the unpo¬ 
larized case. It is straightforward to evaluate the current 
(including the effect of recoil in the phonon-^He scatter¬ 
ing) 


J3 — '^3'^ph dJrec 


dn 3 

dz 


ddx^rec 


dT 

dz 


(24) 


where the ^He diffusion constant, including recoil correc¬ 
tions, is [3 

I^rec — 3^ (i?,C)rVr reci (25) 


and the “thermoelectric” coefficient is 


I^T^rec — {.^1 C) rec- (26) 

In these expressions the basic forms of D and Dt are 
modified by the wall scattering factor 


as shown in Ref. Q the ^He-phonon collision term is di¬ 
agonalized by expanding Wp in Sonine polynomials 





3m* f 
dT 

rec '^3wd 


-1 


( 27 ) 
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FIG. 1: (color online) The diffusion constant, Dree, Eq. (1251) 
as a function of the pipe radius, R, for T = 0.35 K (purple, 
dashed), T = 0.45 K (blue, solid), and T = 0.55 K (red, dot¬ 
ted). We assume here that the ®He wall scattering is diffuse, 

C = 1- 

where 

-1 

C = (28) 

'^3wd Ttujs 

is the fraction of the ^He wall scattering rate that is dif¬ 
fuse. 

To see the effect of scattering of ^He from the walls, 
we take the total wall scattering rate to be simply 


T3wd{R) Reff’ 

where V 3 = ■\/^rJrn* is the mean ^He thermal velocity, 
and Deff = 2i?/3 [see Eq. (IA8I) ] is the effective average 
distance from an interior point to the wall of an infinitely 
long pipe of radius R entering the transport [l3l| . The 
effect of wall scattering on the diffusion constant. Dree, 
for example, is shown in Fig. [T]for C = 1- As we see, the 
effect becomes more important for lower temperatures 
because of the decrease in the phonon density. As i? —>■ 
00 , Dree approaches the result without wall scattering, 
which in the vicinity of the operating temperature regime 
of the nEDM experiment, T ^ 0.45 K, is 

Dree = Cm^/s. (30) 

We also show, in Fig. [51 the corresponding effect of 
phonon-wall scattering on the phonon thermal conduc¬ 
tivity, Eq. 0. 

At temperatures of approximately 0.3K and below, the 
mean free path of a ^He quasiparticle in the bulk medium 
is greater than 1 cm, which is comparable to the radii 
of the pipes considered for the nEDM experiment. In 
this regime, the ^He quasiparticles can lose momentum 
not only by collisions with phonons but also directly to 
the walls of the pipe. This situation is analogous to the 



FIG. 2: (color online) The phonon thermal conductivity, Kph, 
from Eq. 0 for T = 0.35 K (dashed), 0.45 K (solid) and 
0.55 K (dotted) as a function of pipe radius. For R <C Iviaa 
the conductivity rises as R (Eq. (|9ll), while for R 3> Ivisc, it 
rises as R^ (Eq. (fTUll 'l. 

(Knudsen) flow of a low-density gas in response to a pres¬ 
sure gradient, except that in the dilute helium solutions 
the force driving the flow of ^He has two components, 
one due to the ^He pressure gradient and another due to 
the collisions with phonons. When the ^He distribution 
function is stationary, the two contributions to the force 
are equal and opposite. The ^He distribution function 
in principle depends not only on direction in momen¬ 
tum space but also on the radial coordinate in the pipe. 
However, the time scales for smoothing out the radial de¬ 
pendence via diffusion or ballistic transport are of order 
several milliseconds, and therefore short compared with 
the overall evolution time scales of the system. A de¬ 
tailed study of this regime lies outside the scope of the 
present article. 

IV. TIME EVOLUTION OF THE 
TEMPERATURE 

We now estimate, using the heat diffusion equation, 
the timescale for heating the fluid. We assume, as above, 
that the heat is carried by the phonons (see also Eq. (89) 
of Ref. @1 and discussion following) and that the relative 
temperature variation and, hence, the variation of Kph 
is small, so that 

de dT 

— =3Sph-^ = -y-Q = KphV^T, (31) 

where e is phonon energy density. Within a few scattering 
times after the application of heat at one end of the pipe 
(z = 0), the temperature there is approximately fixed at 
To -I- AT. We assume that the temperature at z = L, the 
other end of the pipe, is kept constant at temperature Tq 
by a refrigerator. 

To solve the heat diffusion equation it is sufficient to 
consider the average of the temperature, T{z, t), over the 
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cross-section of the pipe, thus avoiding having to take 
into account details of the counterflows within the pipe. 
The solution is given in terms of the modes in the pipe 
that vanish at 2 = 0 and L, sin/c^z, where = vk/L, 
with u here a positive integer: 


T(z, t) = To + Ar(l — z/L) + ^ Cue sinkuZ. 

(32) 

We denote the thermal diffusivity by 

Dth = Kph/^Sph, (33) 


and recognize ?)Sph = (27r^/15) {T/s)^ as the ^He specific 
heat. The condition that T{z,t = 0) = To except imme¬ 
diately at z = 0, implies that the mode weights are given 
by Cu = —2ATli'Tr. The characteristic time, Tth, to set 
up a steady state phonon wind is essentially that of the 
iz = 1 mode. 


_ _J_ _ L^Kph 

D,hkl ?,^^Spn ^ ^ 

For typical conditions in the experiment, 5 mW of heat 
in a 3 cm diameter, 100 cm long pipe at T = 0.45 K, 
the phonon thermal conductivity is 2.4 x 10® erg/s-cm-K, 
T-th ^ 11 nis, and AT = 3 mK. 


V. TIME EVOLUTION OF THE ®He 
CONCENTRATION 


To begin examining the ®He concentration, we con¬ 
sider its steady-state distribution in the presence of a 
heat current or phonon wind. Because, as we shall see 
below, the term involving Dt is relatively small for low 
concentrations, the condition that the ®He particle cur¬ 
rent, Eq. (I24L vanishes, is simply 


D 


rec 


dns 

dz 


'^ph^3 5 


which has the solution 


nsiz) = = n-i^ao{z), 


(35) 


(36) 


where we define the scale height, h = Drec/vph (for the 
example parameters above. Dree = 225 cm^/s, Vph = 
17 cm/s and h = 13 cm), and 


~ _ no L 

- ^L/h _ 1 ’ 


(37) 


with no the initial uniform ®He density. We note that 
the relative size of the term involving Dt^tcc is simply 
the ratio of the scale heights of the concentration and 
the temperature, 

Dr^rec \dT/dz\ _ \dlnT/dz\ _ DreeSph 
Drec\dn3/dz\ \d\TLno/dz\ Kpu 


about 1/1000 for the example parameters given above. 

In the nEDM experiment, the ®He in the system depo¬ 
larizes in time, primarily due to interactions with walls. 
The depolarized ®He will be removed by a phonon wind 
before the system is recharged with more highly polarized 
®He. As above, we consider the simple situation of a long 
pipe with a heater at z = 0 and closed ends. The evolu¬ 
tion of the ®He is governed by a competition between two 
processes: the phonon wind, which were it to act alone 
would push all the ®He to the downstream (large z) end 
of the pipe, and diffusion of the ®He, limited by scatter¬ 
ing with the phonons, which allows the ®He to drift back 
towards smaller z. 

This evolution of the ®He concentration in the presence 
of a phonon wind is described by the diffusion equation 
resulting from Eq. (1^ . 


dn-i{z,t) dns 

—§^ + ”'”‘- 37 -^’' 


d'^ns 
' dz^ 


= 0 , 


(39) 


where we have dropped the Dt term in Eq. . Once a 
steady phonon wind, with a small temperature gradient, 
is established, we may neglect the temperature depen¬ 
dence of Dree and take Vph and Dree to be constant. For 
a pipe with a large length to diameter ratio, we may treat 
the problem as one dimensional, averaging over its cross 
section as we did above for the heat flow. The boundary 
conditions are that the ®He current, ja, Eq. ( 1211 ) . vanishes 
at the two ends of the pipe. 


dnp _ ^ 


(z = 0,T). 


(40) 


To solve Eq. (I39|) with constant Vph and T), we write 
the ®He density as e^^'^^n{z,t) and decompose n{z,t) as 
a sum of time dependent modes nu{z,t) periodic in 2L 
(see Appendix B): 


no{z,t) = e^/2'‘^h.(z,t), (41) 

where n satisfies the boundary condition 

The spatial parts of the mode functions hu (z, t) are the 
complete orthonormal set 


^z/2h 


(j)u{z) = < 


h (e^!^ - 1 ) 


Tj2- 




cos kyZ -I- 27 ^ sin kuz'^ , 


17 = 0 


17 > 1 


(43) 


with ku = 7rz//T and = [(T/2) (l -|- l/{2hku)‘^^] 
The time dependence of the modes is where 


1 


Tu 


0 , 


17 = 0 

(44) 


klD + vU^D = {kl + l/4/i2) D, 17 >1. 
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FIG. 3: (color online) The ®He concentration, X 3 , from the 
solution of Eq. (1391) as a function of z, the distance along 
the pipe, for various times: t = 0 (dotted), 1 (dash double 
dot), 3 (dash dot), 5 (short dash), 8 (long dash) and 20 s 
(solid). The result is shown for typical parameters in the 
nEDM experiment: = 10“^'^ and 5 mW of heat into a 3 

cm diameter, 100 cm long pipe at a nominal temperature of 
0.45 K. 


The solution of Eq. (1391) for an initially uniform density 
713,0 is then 


7 ^ 3 ( 2 :, i) =713,0 (45) 

1^=0 


where 


f {n{z,0)/nl)(j)^{z)dz 
Jo 


= < 


h - l) 
ihuu 


1^72 ’ 


v = Q 


(46) 


. 1 + (^hky) 


(l + . V> I 


FIG. 4: (color online) The ®He concentration, 2 : 3 , from the 
solution of Eq. dSU as a function of time, t, for various posi¬ 
tions along the pipe: the curves correspond to 2 = 0 (dotted), 
20 (dash double dot), 40 (dash dot), 60 (short dash), 80 (long 
dash) and 100 cm (solid). Note that it takes about 5.5 time 
constants, ri, for the distribution at the hot end of the pipe 
(2 = 0) to reach equilibrium. The result is shown for typi¬ 
cal parameters in the nEDM experiment: 2:3,0 = 10“^°, and 
5 mW of heat into a 3 cm diameter, 100 cm long pipe at a 
nominal temperature of 0.45 K, for which n = 1.8 s. 


h, is substantially smaller than the pipe length. Figure [5] 
plots the difference between X 3 and its steady state value 
for several points along the pipe, showing that, after a 
few seconds, the lowest mode, with ti = 1.8 s, dominates 
the time evolution throughout the pipe. 

The results for the evolution of the ^He concentration 
presented here are equally applicable to a heat flush ex¬ 
periment bein g c arried out at Harvard at natural ^He 
concentration |14| . There one must use the more general 
phonon thermal conductivity as derived in phonon- 
wall scattering in this regime plays a negligible role. 


VI. SUMMARY 


As we show in Appendix B, the general solution may 
be written in compact form in terms of a Green’s function 

fi{z,t)= f g {z,z',t)n{z',0)dz', (47) 

Jo 

where 

00 

g ( 2 , 2 ', t) = ^ ( 2 ) (f)^ ( 2 ') (t). (48) 

The 2 and t dependences of 2:3 are shown in Figs. [3] 
and m respectively, for the case of the uniform initial 
distribution and typical experimental values (5 mW heat 
into a 3 cm diameter, 100 cm long pipe at T = 0.45 K). 
As the figures illustrate, the concentration scale height, 


We have calculated the transport properties of dilute 
mixtures, 2:3 10 “®, of ^He in superfluid "‘He at tempera¬ 

tures around 0.5 K where phonons are the dominant exci¬ 
tations of the superfluid. In this regime, we considered a 
simple one dimensional geometry (a pipe), a heat current 
generates a phonon wind with phonons in local equilib¬ 
rium corresponding to the temperature at that point in 
the pipe. On the other hand, phonon scattering distorts 
the ^He distribution from equilibrium. Starting from the 
known phonon-phonon and phonon-^He scattering, we 
calculate the transport coefficients from the Boltzmann 
equation. We show that, in the presence of a heat current 
which generates a temperature scale height much larger 
than the length of the pipe (i.e., a small relative tempera¬ 
ture gradient), the scale height for the ^He concentration 
can be much less than the pipe length. This leads to 
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tis) 

FIG. 5: (color online) The absolute value of the ®He concen¬ 
tration difference, \x 3 — 2 : 3 , 00 !, from the solution of Eq. (I22II 
as a function of time, t, for various positions along the pipe: 
the curves correspond to z = 0 (dotted), 20 (dash double 
dot), 40 (dash dot), 60 (short dash), 80 (long dash) and 100 
cm (solid). We note that the lowest mode, corresponding to 
the time constant ri, dominates the time evolution after a 
few seconds. The result is shown for typical parameters in 
the nEDM experiment: 2 : 3,0 = 10”^'^, and 5 mW of heat into 
a 3 cm diameter, 100 cm long pipe at a nominal temperature 
of 0.45 K, for which n = 1.8 s. 


a large decrease in concentration at the hot boundary 
and a corresponding increase at the cold end. For tem¬ 
peratures below 0.5 K the mean free paths of both the 
phonons and the ^He can reach the centimeter scale; in 
these cases, scattering from the walls of the container 
becomes important. Finally, we calculate the timescales 
associated with the evolution of both the temperature 
and concentration distributions; because of the large su¬ 
perfluid thermal conductivity, the thermal timescales are 
on the ms scale, whereas the corresponding scale for evo¬ 
lution of the concentration is on the scale of seconds. 
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Appendix A: Phonon drift velocity in the ballistic 
limit 

Here we derive the spatial dependence of the phonon 
drift velocity when the phonon-phonon scattering mean 
free path for large angle scattering, Eq. (O, is much larger 


than the pipe radius, R. For the pipe geometry consid¬ 
ered in the text, this condition holds at a temperature 
r;$,0.3 K. We assume that the pipe axis is in the z di¬ 
rection, and that the transverse coordinates are x and y. 
We also assume that a phonon striking the cylinder wall 
is diffusively reflected, with a distribution of final mo¬ 
menta given by the local temperature, T(z) = Tq + T'z, 
where T' < 0 is the temperature gradient. Then n^{r), 
the number of phonons of momentum q at point r, is 
given by the equilibrium distribution n^{z') at the point 
r' = [x',y',z') on the pipe wall where the phonon at r 
originated: 


r(f) = 


gSg/T(.'(9)) _ 1 - ”9 ^Todq’ ^ ’ 


to lowest order in T', where is the equilibrium distri¬ 
bution. 

The point of origin is determined by simple geometry, 
namely, f— r' = qJ), where T) = |r — r'\. We measure q 
in polar coordinates 9^ and (jjq. Then z' = z — 27 cos 0^, 
X = x'—T) sin0q cos (j)q, and y' = y—T) sin0q sin^g. Using 
xi 2 _|_ yi 2 _ j ^2 cylinder wall, we have then 

sin^ 6q — 2p'Ds\n0q coa{(j)q — cj^r) — = R^, 

(A2) 


where 4 >r is the azimuthal angle of f and p = 3 /x'^ + y^ ■ 
The solution is 


V = 


sin 6 a 


pCOs{4>q - 4>r) + \J- p"^ ain^ {(j)q - 4>r) 


(A3) 


without loss of generality we take = 0. 

The local phonon flow velocity Vph,z[p) is given by the 
local total momentum flux density in the z direction di¬ 
vided by the normal mass density of the phonons, Pp^,, 

1 f d^q 

Vph,z{p) = — 77r^(icoa0qn^{r). (A4) 

Pph J 


To first order in T', only the 27 term in the distribution 
function survives the angular average in the numerator, 
so that 


^ T'J d^q coa^ 6qq'^V driq/dq 

VphAP) = j^3qq2 dnOJQq 


(A5) 


independent of z. Since 27 is independent of q, the inte¬ 
grals over q in numerator and denominator cancel, and 


'>^ph,Ap) 


-is^{V cos^ 9q) 
20 


(A6) 


where the angular brackets denote the average over an¬ 
gles of q. 

The flow velocity averaged over the cross section of the 
pipe (denoted by an overline) is simply 


_ I T' 

'^ph,z ^^2 J P ^P ^ph,z{p') f ■) 

(A7) 












where 


1.0 


1 

= 2tt p dp {V cos^ 9q). (A8) 

In evaluating 'Deff-, the integral over 6g decouples from 
those over p and 4>q, and the latter are easily performed 
if one integrates over p before integrating over tpq- One 
finds 

'Deff = ^ = (A9) 

which expresses the fact that the leng th important for 
thermal conduction is one half of (2?), the average dis¬ 
tance to the wall of the pipe, averaged over the cross sec¬ 
tion of the pipe. Since Vph,z = —i{Kph/TCph)T', where 
Cph is the phonon heat capacity, we find the phonon ther¬ 
mal conductivity, 

Kph = sCphVeff = -sCphR, (AlO) 

which is the Casimir result, Ref. [HI- 
Locally, 

y ^^cos^^g = (All) 

where the elliptic integral. 



must be done numerically. As we see in Fig.[6]the velocity 
profile 

l7rs7?T' 

Vp,Ap) = -(A13) 

is independent of z and nearly quadratic almost to the 
edge of the pipe, where it falls more rapidly, but unlike 
when viscosity dominates, it does not fall to zero at the 
pipe wall. 


Appendix B: Exact solution of the time dependent 
diffusion equation with advection 

Here we construct the general solutions of the time- 
dependent diffusion equation (|3^ by first transforming 
the equation into self-adjoint form by writing 713 ( 2 , t) = 
e^^^^h{z,t). As a result, n{z,t) obeys 

in the interval 0 < z < L (we simply write a generic 
diffusion constant, D, to simplify the notation). The 
boundary condition of vanishing current at the two ends 
of the pipe then becomes 

dh h 
dz 2h 


0.9 


> 0.8 


0.7 


0.0 


0.2 


0.4 0.6 

ip/Rf 


FIG. 6: Normalized velocity distribution across the cylinder 
as a function of (p/R)'^ in the Casimir limit. 


where h = D/vph- To realize the boundary conditions 
we expand h in a complete set of normalized solutions of 
Eq. (IBII) that are periodic in the interval 0 to 2L, 


4>viz) = ctv 1 COS kz -\——- sin kz 
' 2hk 


with k = -Kv/L (y > 1), and 

a,= [{Ll2){l + ll{2hk,f)]-^'\ 

together with the stationary solution 


(B3) 


(B4) 


= 


^z/2h 


[h {e^/h-1)]^/^' 

The modes 4>u{z) form an orthonormal set obeying 


(B5) 


/ (j)y{z)cj)^i{z)dz = 5yy 

Jo 


(B6) 


as well as the completeness relation in the interval 0 < 
2, 2' < L, 


'^(l)^{z)(t>u{z) = 5{z-z ). (B7) 

1^=0 


The time dependence of the modes (j)v{z,t) is e 
with 

and I/tq = 0. Then 


with 


n{z,t) = ^^h^(j)y{z)e 
1^=0 


U - 


77.(2, 0)(j),z{z)dz. 


(B9) 


(2 = 0,L), 


(B2) 


b, 


0 


(BIO) 
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We can thus write the solution n{z,t) in terms of the 
initial density distribution as 

h{z,t) = ( Q{z,z',t)n{z',0)dz', (Bll) 

Jo 

where 

OO 

g{z,z\t) = ^(l)j,{z)(l)j,{z')e~^^^''0{t) (B12) 

is the Green’s function for the diffusion equation in the 
form (EH), with 6 the Helmholtz unit step function. At 
t = 0, the completeness relation implies G{z, z', t) = 5(z— 
z'), and thus 

in the interval 0 < z < L. 


We now convert back to the original form of the diffu¬ 
sion equation, dSH, and have 

n 3 iz,t) = [ G{z,z',t)n3{z',0)dz\ (B14) 

where 

G{z,z\t) = e^^+^''>/^^g{z,z',t) (B15) 

is the Green’s function for the original diffusion equation 

( 15 ^ 1 ) . i.e., 

/ 5 d 5^ \ 

= e(^+^')/2'*<5(z-z')^W- 

(B16) 
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